Dissipation in a rotating frame: master equation, 
effective temperature and Lamb-shift 
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Motivated by recent realizations ol microwave-driven nonlinear resonators in superconducting 
circuits, the impact of environmental degrees of freedom is analyzed as seen from a rotating frame. 
A system plus reservoir model is applied to consistently derive in the weak coupling limit the master 
equation for the reduced density in the moving frame and near the first bifurcation threshold. The 
concept of an effective temperature is introduced to analyze to what extent a detailed balance 
relation exists. Explicit expressions are also found for the Lamb-shift. Results for ohmic baths are 
in agreement with experimental findings, while for structured environments population inversion is 
predicted that may qualitatively explain recent observations. 
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The prospect to tailor quantum devices on ever grow- 
ing scales has stimulated major experimental research in 
the last years. In solid state physics various types of non- 
linear resonators have been fabricated such as mechanical 
beams [l| and superconducting tunnel junctions embed- 
ded in cavities [2f, partly to study fundamental physical 
phenomena, partly to develop highly sensitive amplifiers 
for detection schemes towards the quantum limit. These 
systems can be tuned over broad ranges of parameter 
space and particularly between the domains of classical 
and quantum behavior since they interact inevitably with 
surrounding degrees of freedom. Theory is challenged as 
now dynamical processes like bifurcation, parametric am- 
plification, and period doubling appear on a level where 
h tends to play a decisive role. 

A particular example of this type are microwave-driven 
Josephson junctions (JJ), recently realized in form of the 
Josephson Bifurcation Amplifier [H, Q and the Cavity Bi- 
furcation Amplifier [2]. They are based on the fact that 
near the first bifurcation threshold the switching between 
the two driving-induced states is extremely sensitive to 
parameters of the J J. The latter ones may be influenced 
by the coupling to a quantum system of interest, e.g. a 
Cooper pair island acting as artificial atom [f|, meaning 
that the state of the quantum system can be retrieved 
from the measurement of the switching rate of the JJ. 
It has been shown that this strategy allows for a sin- 
gle shot read-out close to the quantum non-demolition 
(QND) limit @. Theoretically, a general approach is 
provided by the Floquet representation @, but in the 
above situation a more powerful procedure for analyti- 
cal investigations is to describe the driven dynamics in a 
frame rotating with a frequency equal to the response fre- 
quency of the system as already analyzed in the classical 
regime by Dykman and co-workers [H, 0. The exten- 
sion to the quantum regime has been given in [10Hl2j . 
In essence, one arrives at a time-independent description 
with a non-standard Hamiltonian though. A complete 
understanding, however, must include also the bath de- 
grees of freedom residing in the lab frame. 

The common formulation for dissipative systems in the 



weak coup ling regime is provided by so-called master 
equations [13L Hij|, It is well established that the im- 
pact of the bath is twofold, namely, on the one hand to 
induce transitions between eigenstates of the bare system 
so that they are asymptotically populated according to 
the temperature of the bath, and on the other hand to 
shift oscillation frequencies of the off-diagonal elements 
of the density matrix, an effect known as Lamb-shift. 
For driven systems the concept of an effective tempera- 
ture has been introduced in [10j, but it was shown only 
in certain limits [ll|, that this temperature signifi- 
cantly differs from that in the laboratory. In these and 
related works [IH, [l6| the focus has been on the switch- 
ing between driving induced states, while an analysis of 
the structure of the dissipative dynamics in the rotat- 
ing frame has seen less attention. Here we fill this gap 
in consistently deriving a corresponding master equation 
and, importantly for ongoing experimental activities, in 
providing an analysis to what extent a unique effective 
temperature can actually be introduced. This issue is 
intimately related to the existence of a detailed balance 
relation in the rotating frame. 



I. STANDARD FORMULATION 

In the microscopic model for quantum dissipation a 
quantum mechanical system is coupled to a thermal bath 
so that the total Hamiltonian takes the form H = Hs + 
Hb + Hj with Hs being the Hamiltonian of the bare 
system, Hb the Hamiltonian of the heat bath, and Hj 
the interaction, i.e., 
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The time evolution of the density matrix of the full 
compound Wit) obeys the Liouville-von Neumann equa- 
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tion ihdW(t)/dt = [H, W(t)] with an initial state W(0). 
The relevant operator is the reduced density p(t) — 
trs{T / K(t)} for which a simple equation of motion does 
in general not exist. In case of weak friction and suffi- 
ciently fast bath modes, however, progress is made within 
a Born-Markov approximation. One then obtains the 
usual master equation which we cast in the form 



(2) 



where 



C qq [p] = / dsK'(s)[q,[q(-s),p(t)}} 
Jo 

dsiK"{s) [q : {q{-s),p{t)}\ 



(3) 



contains the position operator q(s) in the interaction rep- 
resentation and {,} denotes the anti-commutator. The 
effective impact of the bath appears as a force-force cor- 
relator K{t) = (*)£((% = K' + iK" with £ = 



K(t) 



"dio 
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coth I - | cos(u;i) — z sm(wi) 
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including the spectral density I(u>) = -| J2 n m t ~ 
uj n ) of the bath modes. Processes mediated by these 
modes are revealed most clearly in the representation of 
the master equation using the eigenstates of the bare sys- 
tem Hs\n) = E n \n). An additional rotating wave (or 
secular) approximation (RWA), where off- resonant (fast 
oscillating) terms in © are neglected, gives 
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First, the reservoir induces transitions between eigen- 
states with rates 7„ m = |(?i|<7|to}| 2 D(oj nm ) at frequencies 
<*W = (E n - E m )/h where 
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)e^ s = hl(u)np(u) (6) 



with np(uj) = l/[exp(/37kj) — 1]. The decay rates for the 
off-diagonal elements of the density matrix 
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'-((m\q\m) - (n\q\n)) 2 (7) 

contain the collective rates T mm — 7fc m - Second 

m 

Auj nm by the so-called Lamb-shift 
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oscillation frequencies are renormalized w nm = oj nm 



1/ I l; \,2 DL ( UJ km) , s 

-|M#)I 2H + ( m ^ n > 



(8) 



The bath appears in form of 



D l (uj) = 21 dt K'(t)sw(ut)-K"(t)cos(u)t) (9) 
Jo L - 1 

= 2i^(«) - M[urkta(u) - 77(0)] (10) 



determined by 
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Moreover, M denotes the mass of the system particle and 
?7 S i n is the sin-transform of the classical friction kernel 
r)(t) = J du)I(uj)cos(u)t)/Miruj. 

Since the environment rests in thermal equilibrium its 
probabilities D(u>) to emit and D(—lu) to absorb photons 
are related by detailed balance, i.e., 



oj>0. 



(12) 



This latter condition can also be used to define the tem- 
perature /3 = l/h^T imposed asymptotically (for long 
times) by the bath onto the thermal state of the system: 
the stationary populations P n — p nn of the system obey 
Pn/Pk = P > {— LjJ n k)/D(<jj n i t ), thus leading to the known 
Boltzmann distribution. 



II. MASTER EQUATION IN A ROTATING 
FRAME 



We now consider a system with a standard Hamilto- 
nian Hg = p 2 /2M + V(q) where in 



V(q) 
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Vi(q) 



a harmonic term is splittcd off from anharmonic con- 
tributions. This system is subject to external periodic 
driving with frequency ujd . Then, asymptotically the re- 
duced density does not reach thermal equilibrium, but 
rather approaches a time dependent state with periodic- 
ity Tp = 2-k/lof. For instance, the standard dipole cou- 
pling between system and driving force oc <7COs(iL>i)i)[3] 
leads to ujp = u>d, while a parametric driving of the 
form oc cos(oj£>i)g 2 [Io| gives rise to period doubling 
ujp = Ud/2. As already noted above, particularly the 
former situation is relevant for mesoscopic detectors built 
with microwave-driven J J. The domain where the detec- 
tor is most sensitive, is located around the first bifurca- 
tion threshold, where the anharmonicity V\ is still weak 
and (ujq — lod)/uj{) <C 1. When amplitude and/or fre- 
quency of the external drive are tuned into this regime, 
after a transient period of time the density takes the form 
p(t) ~ p(t) cos(wp-f) with p changing only slowly in time, 
i.e., on time scales t 3> 1/u>f- In principle, to derive 
the master equation for p, one may follow the strategy 
outlined above, however, only after switching to a mov- 
ing frame. We mention in passing that this mapping is 
closely related to a representation within Floquet theory 
7] combined with a projection onto the subspace of the 
leading harmonics ~ exp(±icjpi). 

For this purpose, we introduce a unitary operator of 
the composite system 



U(t) = U s (t)U B (t) = e 



(13) 
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where a and b n are annihilation operators for harmonic 
oscillators in the system and in the bath, respectively. 
The total Hamiltonian in the rotating frame 



H = U f 

follows upon 
exp(±ifccjj?i) , fc| 



U = H S + H B + H I 



discarding fast oscillating terms 
> 1 as a time-independent Hamil- 



tonian of the form Hs(Q, P) — UgHs Us and 
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The operator "Hs = H$ — ihd t coincides with the ex- 
tended Hamiltonian in the Floquet description so that 
the spectrum of H$ reproduces the corresponding Flo- 
quet quasi-energies. Its eigenstates \tp n ) are related to 
the Floquet states \ip n ) via \ip n ) — PfUs^u) with Pp be- 
ing the projector onto the (k- WF)-subspace with k = ±1. 
We emphasize that the original Hamiltonian (p} cannot 
be regained from (|14p in the limit ujp — > since then 
fast oscillating terms neglected in the derivation of (fT4|) 
contribute. 

For the derivation of the master equation the specific 
form of Hs is not relevant. Important are the new bath 
parameters 



1 - UJp/uJn 
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Further, in the rotating frame the system-bath coupling 
is not just a position-position interaction but includes 
also the momenta, a point which we will discuss in more 
detail below. Starting now with ihdW/dt — [H, W] and 
following the standard procedure, one obtains the master 
equation in the rotating frame 



dp 
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Here operators Cqq and Cpp are defined according to 
(|3]) with q replaced by Q, P/(u)fM) with the force-force 
correlator defined by 

K — K' + ? K" 

= \ {F x (t) F y (0)>£ with x, y = Q, P (16) 



where 



C n X Ti 



*p = £c„^-. (17) 



Whereas the mixed operators are 
1 



£qp[p\ 



ujpM J 
1 



ujpM J 



dsK' QP (s) [Q,{P(-s),p(t)}ll8) 
dsiK'^ P (s) [Q,[P(-s),p(t)]\ 



and Cpq is defined according to (fl~8f with Q o P. The 
force-force correlators can be expressed as 



K QQ (t) = K PP {t) 



TT I. 



coth 



(lu + bJp)h/3 



and 



K QP (t) = -k PQ (t) = 



-/(«){ 



coth 



(lu + uip)hf3 



(19) 

cos(wt) — isin(cj£) j 
(20) 

sin(cji) + i cos(cj£) j 



with the spectral density in the rotating frame 

I(uS) = I{lu + lu f )/4. 

Note that the unitary transformation (|13p does not affect 
the equilibrium density of the bath since [Up,, Hp] = 0, 
but just the dynamics of the correlator. Accordingly, 
equilibrium properties in (|20[) appear through a fre- 
quency shift lu — > lu — ujp meaning that in the rotat- 
ing frame the bath carries modes with "negative" fre- 
quencies. The above analysis shows that in the rotat- 
ing frame a typical friction strength is given by i) — 
I(ujp) / AMlof which is smaller than that in the lab frame 
rj = lim w _,.o I(lu)/Mlu. The frequency shift in the bath 
correlation function K(t) effectively produces a decay in 
time which qualitatively is similar to that of K(t). The 
condition for the validity of the master equation (|15[) can 
thus be estimated to read fjh/3 <C 1 provided a typical 
bath cut-off frequency lu c sufficiently exceeds r/ > fj. 

It is instructing to express the above master equa- 
tion in terms of annihilation and creation operators. 
This is most conveniently done in the interaction pic- 
ture representation where the reduced density is given 
by p r (t) = e lflst / h p{t)e~ lilst l h . This way, CE]) trans- 
lates into 



.^dpi 
in—— 

dt 



i C QQ + C QP + C PQ + £pp) \Pl] 



(21) 



We now introduce operators 



aH = J2 \E)(E\a\E')(E'\ (22) 

E'-E=hu> 

ot( w ) = Yl \E)(E\J\E')(E'\, (23) 



E-E'=hui 



Lu n m n 



where the sum runs over all energy eigenstates \E), \E') 
of Hs with fixed energy difference Hlu. a and a' are anni- 
hilation and creation operators of a harmonic oscillator 
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with frequency Q, i.e. 



MQ 
~~2h 

MO 

~2h 



Q 



Q 



iP 
MQ 

iP 



(24) 



The above frequency dependent operators satisfy 
[Hs, a H UJ )) = ftua){ui) and [Hs,a(u>)} — —huja(uj) and 
are thus creation and annihiliation operators at frequency 
u). Note that these relations hold independent of the fre- 
quency Q in the definition (|24[) . For a purely harmonic 
system with frequency ujq, however, the a'(u), a(u) re- 
duce to a, a* only if f2 = cjo- It is thus convenient to 
chose O as the frequency of small oscillations around one 
of the stable extrema of Hs, f° r instance that one where 
the dynamics starts initially. 

The sum of dissipative operators in (|2ip can now be 
expressed as 



QP 



pp 



2MO 



(25) 



with 
L 



& —it(u-u') x 



u,u>' [Pi] 

{k' QQ {u')[a{u>\ [A- a V) + A+a(-u/), pi( 
+ ip J K' QP (u)')[a(cj),{at(u>'),p I (t)}} 
+ iK£ Q (u')[a(uj), {A- ot(a/) + A+ a(-cj'), p z (t)}] 
- nty P {w')[a{u), [aV),p/(t)]]} + ft-c (26) 



Here, the Laplace transforms of the damping kernel ([201 
follow from 



KM 



dtK' xy {t)e 

i 

dtK'Uty 



(27) 



and the impact of the momentum dependent coupling 
terms between system and bath are taken into account 
by coefficients 



2fi , n 2 

/l = . X ± = 1 ± -5- 



(28) 



As for the standard further simplification arises 

within the RWA, where all terms in the sum ([26)) with 
U =/= uj' are assumed to oscillate so rapidly on the time 
scale for relaxation that they may be discarded. Accord- 
ingly, (|26l) reduces to 



£oo + £ P q + £qp + C PP | RWA - ^-g^ X! Lu ' u 



We recall that the anharmonicity of the system po- 
tential gives rise to the frequency sum. In contrast, 
for a purely harmonic system with frequency 57 = d>o 
it collapses to contributions for cj = ±Q. Moreover, 
a(— = a'(— \u>\) = so that terms carrying in ([26]) 
the coefficient A + vanish. One should emphasize that this 
is not the case if in (|2~4"]) a frequency fi ^ uiq is chosen. A 
conventional master equation then emerges with renor- 
malized friction functions K' — > (\~K'qq — [iK'qp) and 

K" — > (X~Kqq + iiK'qp). In the Wigner representation 
the momentum coupling terms ensure that in the labora- 
tory frame the known Fokker-Planck equation for driven 
systems is reproduced. To neglect them in the rotating 
frame leads in turn to spurious diffusion terms in the lab 
frame [l7| . 

In the anharmonic case and for elevated temperatures 
it may be more convenient not to work explicitly with 
annihilation/creation operators, but rather use a repre- 
sentation as in (0, where matrix elements of those system 
operators appear which couple to the bath. In particular, 
these matrix elements can then directly be evaluated for 
the full anharmonic problem by means of semiclassical 
techniques in certain ranges of parameter space [ToL [l8| . 
The standard result can easily be generalized to the ro- 
tating frame situation: in addition to matrix elements 
(m\Q\n) also (m\P/Mcjp\ri) must be taken into account; 
these matrix elements are multiplied by the bath corre- 
lation functions D xy and D xy derived from K xy . The 
diagonal part of the density p thus obeys a Pauli-master 
equation with properly modified transition rates, 



:\(n\P\m)\'< 



Inm = D QQ (u} mn ) \(n\Q\m)\ 



(uj f M)' 



DQp(uJ m n) 

uj f M 



((n\Q\m)(m\P\n) - (n\P\m)(m\Q\n)){30) 



(29) 



If typical transition frequencies in Hs are small com- 
pared to the driving frequency wp the dominant process 
in the transition rates is given by the position matrix 
element. For instance, close to a harmonic minimum 
with frequency Coq and near the bifurcation threshold the 
momentum dependent terms provide contributions which 
are suppressed by factors on the order of wq/wf < 1. 

As already mentioned above, the RWA does not always 
apply. For instance, for higher lying states the energy 
spectrum of Hs may have accumulation points with a 
dense distribution of transition frequencies. This issue 
becomes particular relevant for the situation described 
above, namely, a driven system close to its first bifur- 
cation threshold. Then typically Hs exhibits two stable 
domains in phase space which are separated by an un- 
stable one, as e.g. for parametric driving where one has a 
double well structure with the two minima corresponding 
to the stable extrema and the barrier top to an unsta- 
ble saddle point of the dynamics. Accordingly, for suffi- 
ciently deep wells the RWA applies for low lying states. 
By tuning the system close to a bifurcation point, how- 
ever, extrema tend to coalesce so that the spectrum of 
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Hs becomes narrowly spaced and the RWA is no longer 
applicable. 



III. EFFECTIVE TEMPERATURE 



The model described by (|T4|) can be considered as a 
system subject to relaxation as in the standard situa- 
tion. Effectively for long times the bath dictates its tem- 
perature to the system such that the populations of the 
quasi-energy levels become stationary and are distributed 
according to a balance between emission and absorption 
processes. In case of a purely harmonic system the cor- 
responding stationary populations can be calculated ex- 
plicitly [121 • F° r the anharmonic case they follow within 
the RWA from the extended Pauli-master equation. It 
is thus convenient to introduce the concept of an ef- 
fective temperature, which is defined according to (fT2")) 
by exp(-/3 xy ?ko) = D xy (-u>)/ D yx (uj), with x,y = Q,P. 
Now, a quantum of energy fko > emitted from the sys- 
tem reaches the bath in the lab frame either with energy 
fiw + = Ti{ujf +lo) or tko^ = H(luf — w), thus determining 
two distinct regimes. In the range to > top only lo + > 
so that only one channel of bath modes is accessible and 
one unique effective temperature follows from the stan- 
dard expression with a shifted frequency 



/8(l + 



UJp\ 
UJ I 




FIG. 1: /3qq (thick) and Pqp (thin) vs. frequency for fixed 
lab temperature fthwF — 4 for an ohmic bath 



(31) 




The situation is different in the second range uj < ujp 
where both frequencies lj± > and two channels in the 
bath are open to give different effective temperatures for 
different diffusion processes, namely, 



QQ 



1 ln [ MfaJ+) + 1] + n p (w-)I(u.) 
flu \ np{u) + )I(w + ) + [np(w-) + 1]I(cj_) 



QP 



1 

fvjj 



ln 



—np(u+)I(u)+) + [np(u-) + l]J(w-) 



32) 



•33) 



and from O and fl2TJl also /3^ Q 



/3< p and 0< 



QP — ^PQ- 

These results are illustrated in ( fig 032]) • Apparently, the 
different combinations of the two channels lead to (3qq ^ 
(3q P so that a unique effective temperature cannot be 
specified, which in turn means that a detailed balance 
relation does in general not exist in the rotating frame. 
Only if the temperature tends to zero, i.e. LOph/3 ^> 1, do 
the individual effective temperatures merge (see figlU so 
the detailed balance is reestablished. 

Let us now discuss the situation in more detail. First, 
while one recovers for vanishing rotating frame frequency 
T xy = T xy = T, for finite cop a discontinuity occurs 
T< y (uj F ) > T> y {uj F ) = T/2. Second, for T = one 
has T xy = 0, while the expressions (|32|) and (|33| predict 
that the common effective temperature Tqq = Tq P — 
(fajj/kft)/ ln[/(a;+)//(a;_)] can be finite. Third, for small 
frequencies u <C lof more transparent expressions can be 



FIG. 2: Pqq (thick) and /3 QP (thin) 
cj = 0.2u_f for an ohmic bath. 



vs. lab temperature at 



derived, i.e., 
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HI(u>f) 

2k B I'(u F ) 



• coth 



hut 



(34) 



T. 



I{LO F )(e huJF0 - l) 2 



QP 



2 k B e^Pfvu} F pI{u F ) + r(oJ F )(e 2huJl ' 13 - 1) 



(35) 

with I' = dl/diu. For an ohmic spectral density / = 
M'yui with I(u>p)/I'(u>f) = this is in agreement 
with experimental observations 19], where one always 
has T o < > T and for T = arrives at T< Q = T< p = 
hup/SkE [11]. Accordingly, the rotating frame system 
behaves quantum mechanically only if fko > k B T xy , i.e. 
u)p ^> lo > LOp/2- In contrast, in the high frequency 
domain the temperature is always reduced. 

For a structured environment the two-channcl- 
processes in (|32|) can lead to counter-intuitive phenomena 
as seen in figs. [3l |U There, we consider the coupling to 
a damped harmonic oscillator (frequency uj p , damping 
strength 7), i.e., 



Mjoj 



{uj 2 



,2 ) 2 + 7 2 w 



2,, ,2 



(36) 



0.0 



0.5 

Ld/iOp 



1.0 



FIG. 3: /3qq (thick) and /3gp (thin) vs. frequency for fixed lab 
temperature fthwF — 4 for a damped harmonic oscillator bath 
with 7/ixjf = 0.5 at uj p = .95ojf (dashed) and ui p = 1.35cjf 
(solid). 4 




FIG. 4: Effective inverse temperature vs. lab temperature at 
uj — Q.2uif and for spectral bath densities as in fig. [3] 



for which I'(u)p) < if top > u p . Hence, ()34ll35j) predicts 
a negative effective temperature physically correspond- 
ing to the fact that absorption becomes more probable 
than emission and a population inversion is induced. Re- 
cent experimental observations of enhanced relaxation in 
a quantronium circuit coupled to a cavity bifurcation am- 
plifier 20] may be explained qualitatively by this phe- 
nomenon. We mention in passing that the Markov ap- 
proximation employed in the derivation of (|15j) remains 
valid as long as the frequency aj p exceeds the typical fre- 
quency for the bare system dynamics SI. 



IV. LAMB-SHIFT 



As discussed above, the reservoir also affects the fre- 
quencies of the off-diagonal elements of the density ma- 
trix. This Lamb shift can now easily be obtained from 




-0.2 



FIG. 5: Bath factor for the Lamb-shift in the lab D L (dashed) 
and the rotating frame Dqq (thick) and Dqp (thin) for 
hf3cj F = 4. 



as 



{E E 

r=0 X,Y£{Q,P/Muj f } 



(m\X\r) 



x (r\Y\m) 
— (m — > n)\ 



h 2 



2ft 2 

n\Q\n)(m\P/(upM)\m) 
(37) 



with the friction function 



Dqq — Dpp 



°dSl 



2K' sin {u>) + 2 —(u F -SI) d (Sl) (38) 

JO 7T 



and 



D% p (lj) = -D PQ (u) 



2 / dt 
'o 



K' QP (t) cos(wt) + K'^ P {t) sin(art) (39) 

id ' J 2hp7T 



n— — oo 

r°°dSl 
+2 / — ujd {Sl) 



o 



7T 



with 



~, . . ^ r°° dSl u 



(SI) 



and 



d n (Sl) 



I{Sl)Sl 2 



SI 2 + vl cu 2 -{Sl- cup) 



(40) 



(41) 



(42) 



and the Matsubara frequencies v n — 2iin/hj3. Since 
for a purely ohmic spectral density the above expression 
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diverges in the low as well as in the high frequency range, 
we take /(cj) = M^uj^uj 2 /(uj 3 +uj 3 ) (seengjSJ). Deviations 
are small in the high frequency domain, but pronounced 
effects occur for w w uip such that even the sign of the 
respective bath factors changes. 

V. CONCLUSIONS 

To summarize we have provided a consistent deriva- 
tion of a master equation needed to capture the dis- 
sipative dynamics of periodically driven nonlinear res- 
onators near the first bifurcation. Various expressions 
for the corresponding master equation have been derived 
and discussed. It turns out that a position-position in- 
teraction between system and bath in the laboratory 
frame translates into additional momentum-momentum 
and momentum-position couplings in the rotating frame. 
In case that typical transition frequencies in the rotat- 
ing frame are sufficiently small compared to the external 
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driving frequency, these latter terms are small compared 
to the dominant position-position coupling. The concept 
of an effective temperature has been introduced to an- 
alyze to what extent a detailed balance relation exists 
in the rotating frame. In the strict sense, one recovers 
detailed balance only at very low temperatures. For a 
structured environment phenomena such as negative ef- 
fective temperatures are predicted. Explicit expressions 
have also been given for the Lamb-shift, which allow to 
better understand recent and future experimental find- 
ings. 
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